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ABSTRACT 

Context. The chromospheric layer observable with the He I 10830 A triplet is strongly warped. The analysis of the magnetic mor- 
phology of this layer therefore requires a reliable technique to determine the height at which the He I absorption takes place. 
Aims. The He I absorption signature connecting two pores of opposite polarity in an emerging flux region is investigated. This signa- 
ture is suggestive of a loop system connecting the two pores. We aim to show that limits can be set on the height of this chromospheric 
loop system. 

Methods. The increasing anisotropy in the illumination of a thin, magnetic structure intensifies the linear polarization signal observed 
in the He I triplet with height. This signal is altered by the Hanle effect. We apply an inversion technique incorporating the joint action 
of the Hanle and Zeeman effects, with the absorption layer height being one of the free parameters. 

Results. The observed linear polarization signal can be explained only if the loop apex is higher than a;5 Mm. Best agreement with 
the observations is achieved for a height of 6.3 Mm. 

Conclusions. The strength of the linear polarization signal in the loop apex is inconsistent with the assumption of a He I absorption 
layer at a constant height level. The determined height supports the earlier conclusion that dark He 10830 A filaments in emerging 
flux regions trace emerging loops. 
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1. Introduction 

The magnetic field is commonly measured in the solar photo- 
sphere, while measurem ents of the chromos pheric and coronal 
field are rarer (see, e.g.. ISolanki et alj|2006l for an overview). 
Whereas photospheric magnetic field measurements dominantly 
employ the Zeeman effect, in higher layers of the solar atmo- 
spher e other methods come increasingly into play (see iLaggl 
I2OO5L for a review), in particular also the Hanle effect. One of 
the most used diagnostics of the magnetic field in the upper 
chromosphere is the 10830 A He 1 multiplet. For this multi- 
plet the Zeeman effect is sensitive to magnetic fields above ~ 50 
Gauss, well suited to study magnetic field s of active regions (e.g . 
iHarvev & HaUl I197U iRiiedi et alJ 119951: iPenn & KuhnI Il995h . 
The Hanle effect is sensitive to weaker fields (1-100 Gauss) and 
has been used to measure the mag netic fields of prominences 
and spicules observe d off-limb (e.g. iTruiillo Bueno et a 
iMerenda et al.ll2066l) an d of filaments and enierging flux regions 
obser ved on-disk (e.g. ITruiillo Bueno et all et alj 

I2OO4I) . Unlike the level of polarization produced by the Zeeman 
effect, the polarization produced by scattering and then modified 
by the Hanle effect (in both 90° and forward scattering) depends 
on the anisotropy of the radiation field exciting the scattering 
atoms. 

This anisotropy increases with the height of the observe d 
structure above the solar surface jSahal-Brechot et all Il977h . 
Therefore, this effect can in principle be used to determine 
at what height in the chromosphere or corona the observed 
magnetic structures are located. Since the He 1 lines are op- 
tically thin in some locations, notably in the quite Sun (e.g. 



iGiovanelli & Halllll977h . bu t optically thick at others, in partic- 
ular in young active regions (IXu et al.ll2009l) . its formation could 
vary considerably, making such an independent determination of 
the height an interesting addition to the diagnostic capabilities of 
the He multiplet. 
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ISolanki et alJ (l2003h proposed an interpretation of spec- 
tropolarimetric observations of an emerging flux region based 
on the assumption that the He 10830 A line is formed along the 
freshly emerged ma gnetic field li nes, which in this case reached 
a: 10 Mm. Recentlv. lJudgd (l2009h has proposed a different inter- 
pretation of these observations according to which the magnetic 
field sampled by the observations is located in a thin layer at a 
constant height (^ 2.4 Mm) near the top of the chromosphere. 
This controversy lends an even stronger need for an indepen- 
dent diagnostic of the height of the magnetic features sampled 
by the He 1 10830 triplet. In this paper we develop the height 
diagnostic further, first considering the influence of the height 
on the strength of the Stokes Q and U profiles. Then, exploiting 
the sensitivity of the Hanle effect to the height of the observed 
structure, we have analyzed again the data of May 13 2003, in 
order to independently estimate the height to wh ich these data 
refer t o, particularly at the location interpreted bv ISolanki et al.l 
(I2003h to be the apex of the highest reconstructed loop. 

In Sect. |2] we summarize the data we analyze in this pa- 
per and introduce the codes we use to synthesize and invert the 
Stokes profiles. In Sect.[3]we present results of a parameter study 
of the strength of the linear polarization signal. We also compare 
the observed with the theoretical profiles in order to fix a mini- 
mum height and perform inversions with HeLIx"*", which allow 
the field strength, direction and the height to be determined. 
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2. Observational data and codes 

2.1. Data 

Spectropolarimetric maps of the emerging flux region NOAA 
945 1 were recorded on 13 May 2001 using th e Tenerife 
Infrared Polarimeter (TIP. iMartmez Pillet eFani 19991) . The cov- 
ered wavelength range contains the 10830 A multiplet. The ac- 
tive region was located near disk center at yu = cos Q = 0.8. These 
observations have been analy zed in several papers (ISolanki et al.l 
l200llll^ et al.l2004l 120071) which emphasize different aspects 
of this rich data set. In the present paper we focus on one spatial 
pixel of the observed map situated in the central part of the active 
region (pixel jc = 72,y = 43 of Fig. 1 in iLagg et ai]|2004l) . 

The Stokes profi les at this pixel are also displayed in Fig. 3 
of iLagg et al.l (l2004l) . This point is within the area of few pix- 
els where the Hanle effect in forward scattering is most promi- 
nent in the Stokes profiles. Around this area the linear polariza- 
tion rapidly falls below the noise level. This pixel also displays 
the largest linear polarization signal in the entire map and, im- 
portantly, it corresponds to the a pex of the highest loop recon- 
structed by ISolanki et al.l (l2003h . We refer to this pixel as the 
loop apex. It makes sense to concentrate on this loop-top pixel 
since it displays th e largest discrepancy between the interpreta- 
tions of lSolanki et a l. (2003) and Judge (2009) and the large lin- 
ear polarization alsomakes this pixel the best location to demon- 
strate the new technique. 

The degree of linear polarization reaches 0.45% of the con- 
tinuum intensity in the center of the red component of the mul- 
tiplet (see Fig. 3 of lLagg et al.ll2004l) . The observed polarization 
is approximately parallel to the dark Unes visible in the core of 
the red He I component which connect the footpoints of the re- 
constructed loop spanning the emerging flu x region. These elon- 
gated absorption structures (see Fig. 1 of ISolanki et al.l [20031: 
iLagg et ani2004l) are reminiscent of Ha arch filament systems. 
The Stokes V signal in this region is comparable to the noise 
level of 10"^ of the continuum intensity. The surrounding pixels 
show similar Stokes spectra with slightly lower linear polariza- 
tion signal, but enhanced circular polarization signal. 

2.2. Synthesis code 

The Hanle profiles observed on disk are produced by the pres- 
ence of atomic polarization (i.e. population imbalances and 
quantum coherences between the magnetic sublevels pertaining 
to the upper and lower term of the 10830 A multiplet; i.e. 
and 1 respectively), generated proportionally to the degree of 
anisotropy of the unpolarized radiation field coming from the 
underlying layers of the solar atmosphere. This atomic polariza- 
tion, modified in the presenc e of a magnetic field by th e Hanle 
effect in forward scattering dTruiillo Bueno et al]|2002h . selec- 
tively absorbs and re-emits the unpolarized light coming from 
solar disk center, acting as sinks and sources of linear polariza- 
tion. 

reproduce these observed Hanle profiles we 
synthesis code with the s ame forward calcu- 
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lAsensio Ramos et al., (2008^ 7 

The Helium absorption is modeled to take place in a 
constant-property slab located at a given height h above the solar 
surface and illuminated from the underlying solar photosphere. 
In this situation the degree of atomic polarization of the Helium 
atoms is a function of (1) the anisotropic illumination that pro- 



duces this population imbalance, and (2) the presence of a mag- 
netic field that modifies it. In a plane parallel atmosphere the 
anisotropy of the radiation at a given A depends basically on the 
height of the scattering atom. 

In order to quantify the atomic polarization we solve the 
statistical equilibrium equations for the multipole components 
Pq{J, J') of the atomic density matrix for a multi-term atom. 
We take into acc ount coherences between different /-l evels 
(see chapter 7 in iLandi DegFInnocenti & Landolfil |2004|) . the 
anisotropy of the radiation field for a given height (computed 
for the center-to-limb variation tabulated in [Pierce 2000) and 
the presence of a magnetic field. Note that the sophistication 
with which the polarization produced by the Hanle effect is 
computed here is co nsiderably greater than in the analysis of 
ISolanki et al.l ( l2003h . From the calculated density matrix ele- 
ments it is possible to compute the emission coefficients (e/, ex 
with X - Q,U, V) and the absorption coefficients (77/, rjx) for 
polarized light. To calculate the emergent Stokes parameters we 
use the solution of t he radiative equations for a s lab of constant 
properties (see, e.g., lAsensio Ramos et"al ] l2008h taken to lie in 
the upper chromosphere 
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where /q and Xq are the boundary conditions representing the 
radiation entering the slab from below. We tak e /p to be the pho- 
tospheric continuum at yu = 0.8 as tabulated bv lPiercel(l2000l) and 
X() - 0. This solution differs from the Milne-Eddington solution 
employed in earlier analysis of the observed Stokes profiles con- 
sidered here. The slab geometry has the advantage that it allows 
a unique height to be assigned to the center of the absorption re- 
gion. It may also be closer to the case of He I formation in an 
emerging flux region. 

2.3. Inversion code HeLIx* 

The forward model described in the previous section is imple- 
mented in the inversion code HeLIx"*" (iLagg et al.ll2010l) . The 
free parameters of the constant-property slab model are the 
magnetic field strength, inclination and azimuth, B, y and x 
(in the solar reference frame), the line of sight velocity, vlos ^ 
the damping parameter a, the thermal broadening of the line, 
VDopp, the optical thickness of the slab, d, and the height of 
the slab above the s olar surface, h. Following the analysis of 
ISolanki et al.l ( l2003l) we assu me the filling factor to be unity. 
The Pikaia genetic algorithm dCharbonneaull 19951) was used to 
maximize the fitness function /, defined as the reciprocal of 
the squared difference between the measured and the observed 
Stokes vector P {P - I, Q, U, or V) over a wavelength window 
of (-1.8, +0.9) A centered at the red lines of the He I triplet, di- 
vided by the strength s of the observed profiles (sj - Yj - 1 1 

A 

and 5q,[/,v = 2 KG' U, V)(A)\) and the number of free parameters 
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The weighting functions wp were chosen to be unity within 
+0.6 A around the core of the red component of the He I triplet, 
and was set to 0.6 outside this range. To take into account the 
reduction of the intensity level caused by the wings of the pho- 
tospheric Si I at 10827 A we fit this line with a Voigt function. 
The parameters of this Voigt function are determined from an 
average of 100 Stokes / profiles centered around the loop apex. 

In extensive tests we found the Pikaia algorithm to be the 
most reliable technique to determine the global minimum within 
the considered parameter space independently of the choice of 
the initial guess values. Additionally, the random walk conver- 
gence of the Pikaia algorithm allows for a comprehensive error 
analysis in a statistical sense. The fluctuation of the retrieved 
parameters for individual, independent inversions of the same 
Stokes vector naturally delivers a measure for the reliability and 
robustness of the retrieved parameters. 

3. Results 

We have performed computations with both codes introduced in 
Sect. I2.2l and l2.3l These computations have different aims. In a 
first step, described in Sect. 13.11 we obtain an idea of the depen- 
dence of the linear polarization in the He I triplet on the com- 
ponents of the magnetic vector and the height of the slab. Then, 
in Sect. 13.21 we test whether this technique can be employed to 
set lower or upper limits on the height of the scattering atoms by 
comparing with a given observed profile. Finally, in Sect. 13. 31 we 
obtain a best estimate of the height of the sampled field using an 
inversion technique. 

3.1. Dependence of the computed polarization on lieiglit and 
magnetic field 

The way in which the presence of a magnetic field modifies the 
emitted polarization signal depends on its direction and on the 
position of the observed point on the solar disk. For example, 
off'-limb linear polarization emission is maximum in the absence 
of a magnetic field, or for a vertical magnetic field (with respect 
to the solar surface). Its value decreases for increasingly hori- 
zontal magnetic fields. For observations at solar disk center, on 
the contrary, the He atoms do not emit linear polarization in the 
absence of a magnetic field or for a vertical magnetic field. Here 
the intensity of linearly polarized radiation increases for increas- 
ingly horizontal magnetic fields. For observations at other loca- 
tions of the solar disk a behavior somewhere in between these 
two extremes is present. Figure [T] displays the results of an ex- 
tensive parameter survey, which aims to determine the depen- 
dence of properties of the linear polarization on the parameters 
of the magnetic vector and the height of the slab. The computed 
examples refer to yu = 0.8, which corresponds to the location 
of NOAA 945 1 at the time it was observed in H e I. E ach dot 

represents the degree of linear polarization p - ^^^^^ plotted 
versus the polarization angle a = 0.5 tan"'(C//0 (both at the A 
of maximum p) for profiles computed for all possible orienta- 
tions of the magnetic field. The inclination angle with respect to 
the local solar vertical, Ob, is varied from 0° to 180° in steps of 
2° (Ob = 0° for a magnetic field vector pointing outward from 
the Sun) and the azimuth angle in the plane parallel to the solar 
surface, is varied from -90° to 90° in steps of 4° (xb = 0° for 
a magnetic field vector lying in the plane of the local solar verti- 
cal and the observer direction). At every distance above the solar 
surface there is a maximum polarization, corresponding to "ver- 
tical" magnetic fields (less inclined than Van Vleck's angle 54.7° 



to the solar surface normal, red points in the figure). This max- 
imum increases with height. The same increase occurs for the 
"horizontal" field (more inclined than Van Vleck's angle, black 
points in the figure), although in this case the maximum reached 
is lower. Note that irrespective of the height there is no minimum 
polarization. This behavior implies that for a given measured lin- 
ear polarization magnitude a lower limit on the height of scatter- 
ing atoms can be placed, but no upper limit (just based on the 
degree of linear polarization). Furthermore, this lower limit is 
higher for a horizontal magnetic field than for a vertical field. 

3.2. Comparing observed with synthetic profiles for different 
heights 

In Fig. [2] we comp are the linear polarization profile of Fig. 3 of 
iLagget al.1 (l2004l) . pixel x= 72, y= 43 in their Fig. 1, (points) 
with the computed profile displaying the maximum linear polar- 
ization for a set of heights between 1 and 10 Mm above the so- 
lar surface. Each solid profile is plotted versus wavelength and 
corresponds to the height given on the horizontal axis. Each syn- 
thetic profile is the one showing the largest U/Ic among all pro- 
files in the "horizontal" range, i.e., more inclined to the vertical 
than 54.7°, formed at that height (for both synthetic and observed 
profile the reference direction for positive Stokes Q form an an- 
gle of 12° with the nearest solar limb, resulting in Q/I^ ~ 0). 
The requirement that the magnetic field is in the "horizontal" 
range is driven by the fact that the Stokes V signal is comparable 
to the noise level. Additionally, the Stokes V spectra along the 
loop structure change polari ty close to the loop a pex. It is also 
consistent with the finding of lSolanki et al. I (l200l that this loca- 
tion corresponds to the top of a loop. Clearly, for heights lower 
than 5 Mm the anisotropy of the radiation field that excites the 
Helium atoms is not able to produce sufficient atomic polariza- 
tion to reproduce the observed U profile. In particular, a slab at 
a constant height of 2.4 Mm, as proposed bv lJudgd (l2009h . does 
not provide a satisfactory description, since even the strongest 
U profile from that slab fails to achieve the polarization level of 
the observed profile. It falls well short of the polarization in the 
5 most strongly polarized points of the observed profile. 

3.3. Results of inversions 

We carried out a series of inversions (see Section 123] for a de- 
scription of the code and setup) in order to determine the at- 
mospheric parameters whic h best descr i be the observed Stokes 
profiles shown in Fig. 3 of iLagg et al.l (l2004l) . The minimiza- 
tion technique used for the inversions with HeLIx"*" involves the 
genetic algorithm PIKAIA (Charbonneau 1995). This algorithm 
has the advantage of finding the global minimum (i.e. the syn- 
thetic Stokes profile best matching the observed one) within the 
considered parameter space completely independent of any ini- 
tial guess values. However, the path of convergence to this global 
minimum is random and may require an infinite number of it- 
erations. To overcome this problem we repeated the inversions 
100 times, each with a sufficiently high number of iterations 
(n = 500). The result of the inversion run with the best fitness 
out of the 100 independent runs was considered to be the solu- 
tion closest to the global minimum. 

The allowed range for the values of the the atmospheric pa- 
rameters for the individual inversion runs is presented in Tab. 
[T] These ranges represent a carefully selected compromise be- 
tween, on the one hand, stability and speed of the inversions and, 
on the other hand, the full coverage of expected values for the 
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Fig. 1. Polarization diagrams computed for different heights above the solar surface (from 1 Mm in the leftmost frame to 9 Mm in the 

rightmost frame). In every diagram, the y-axis gives the degree of linear polarization p = V2" + U^/I^ the x-axis the polarization 
angle a = 0.5 tan^\U/Q). p represents the maximum in the wavelength range covered by the He I multiplet. The plotted a has 
been sampled at the same wavelength. Red points correspond to "vertical" magnetic fields (Ob < 54.7° or Ob > 125.3°), while black 
points correspond to "horizontal" magnetic fields (54.7° < 9b < 125.3°). 



atmospheric parameters at the loop top. We carried out 100 in- 
dependent inversion runs on the loop top pixel for each fixed He 
slab height ranging from 0.7 Mm to 17.5 Mm. The best fitness of 
the 100 inversions at each height is plotted as a function of height 
in Fig. [3] The figure clearly shows that a global maximum of the 
fitness is achieved for a height of the He slab at the loop top 
of »6.3 Mm. Note that a higher fitness corresponds to a better 
fit. The smoothness of the maximum fitness curve in Fig.[3]indi- 
cates that even relatively small fitness difference s are significa nt. 
The fitness at the height of 2.4 Mm, proposed bv lJudgd (l2009l) is 
marked in Fig.[3]by a diamond. It is significantly lower than the 
best fitness. Additionally we marke d the fitness corresp onding 
to height of 10.5 Mm as retrieved bv lSolanki et~an (l2003h . 

Fig. m shows the fit obtained from the best inversion run for 
a slab centered on a height of 6.3 Mm. Lower height values fail 
to produce a sufficiently strong linear polarization signal. If the 
absorbing Helium slab is located higher than 6.3 Mm, the lin- 
ear polarization signal resulting from scattering would be larger 
than the observed one. The inversion procedure compensates 
these strong signal by increasing the magnetic field strength, 
since a stronger magnetic field leads to depolarization of the 
scattering polarization. However, with increasing magnetic field 



strength the Zeeman signal increases, resulting in a mismatch 
for the Stokes V signal and consequently in a reduction of the 
fitness. The atmospheric parameters from the best inversion run 
are displayed in the fourth column of Tab. [1] The magnetic field 
strength at the loop top is 3 1 1 G, the field is nearly parallel to the 
solar surface (inclination angle y=95°), and the azimuth angle 
confirms the orientation of the field along t he visible structures i n 
the He equivalent width image in Fig. 1 o f lSolankietan(l2003h . 
The up-flow of 1.3 km/s at the loop top confirms the rising flux 
tube scenario presented in that paper. 

The simpler analysis of lSolanki et al.l (l2003h provided a field 
strength that dropped to around 50 G at the loop apex, while 
here we obtain 320 G. T he higher value lies consid erably closer 
to the field strength that IWiegelmann et al.l (l2005h retrieved at 
the apex of the loops closest to the one whose apex we study 
in this paper using a force-free extrapolation from the photo- 
spheric magnetic field map (305 - 615 G). This illustrates the 
importance of using the more realistic treatment o f the Hanle ef- 
fect em plo yed here cornpared to that employed bv lSolanki et al.l 
(I2003h and lLagg et all (l2004 . 

We also have inverted the few pixels showing linear polariza- 
tion located near this one, finding comparable values of magnetic 
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Fig. 2. Observed Stokes Ullc profile (points) of Fig. 3 of iLagg et al.l ( |2004|) compared with the maximum Ullc profiles computed 
for each height (solid line). At each height the observed and computed profile is plotted vs. A. The same observed profile is repeated 
for comparison at each height. The computed Stokes / profile (not shown here) fits the observed intensity. The reference direction 
for positive Stokes Q forms an angle 12° with the nearest solar limb. This reference direction gives Qlh ~ 0. The horizontal bar 
indicates the wavelength range, 10827.9 - 10831.6 A of the leftmost profile. 



field and slightly lower heights, consistently with the fact that the 
pixel discussed here in detail was located at the loop apex. 

Table 1. Ranges for the atmospheric parameters. The fourth col- 
umn shows the atmospheric parameters for the best fit. 



parameter 


min. 


max. 


best fit 




value 


value 


(h = 6.3 Mm) 


S[G] 


5 


1000 


311 


rn 


60 


120 


95 


xn 


-90 


90 


-86 


vlos [km/s] 


-10 


10 


-1.3 


a 


0.05 


0.50 


0.22 


VDopp [km/s] 


5 


25 


6.80 


d 


0.20 


0.90 


0.47 



4. Discussion and conclusions 

We have shown that an analysis of the linear polarization of the 
He I 10830 A triplet using an atomic polarization interpreta- 
tion can be used to estimate the height at which the triplet is 



formed and hence the height at which it samples the local mag- 
netic and velocity field. Lower limits on this height are more 
straightforward to obtain than upper limits. Because the lower 
limit is tighter for larger polarizations, this technique works best 
for profiles displaying large-amplitude Hanle-like Q or U sig- 
nals. 

ISolanki et alj ( l2003h have proposed a reconstruction tech- 
nique to obtain the 3D magnetic field structure of an emerging 
flux region observed near the disk center This technique is based 
on the assumption that the He 10830 A line is formed along 
the magnetic flux tube as it emerges through the chromosphere 
into the corona. The reconstruction was constrained by requiring 
the physical quantities in the loops to be smoothly varying and 
for both footpoints of the loops to lie at the same height. They 
found that the tallest reconstructed loop reached in its central 
part a height of ^ lOMm over the solar surface. Recentlv lJudgd 
(l2009h has proposed an alternative interpretation of this obser- 
vation. He argues that the observed magnetic fields are located 
in a thin layer at a constant height (a; 2.4 Mm) at the top of the 
basically plane-parallel chromosphere. 

In the present pap er we have ana l yzed t he most critical pixel 
in the observation of ISolanki et al.l (l2003l) . namely that at the 
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Fig. 3. Best fitness out of 100 individual inversion runs for fixed 
height values from 0.7 Mm to 17.5 Mm. The optimum fit to the 
observed Stokes vector is achieved for a height of the He slab of 
^6.3 Mm. The diamond and the rectangular symbol indicate the 
fitness at a height of 2.4 Mm and 10.5 Mm. 



apex of the highest loop, with the help of HeLIx"*", an inversion 
code incorporating both Hanle and Zeeman effects. We found 
the best inversion fit for an hei ght of =a 6.3 Mm. Th is value 
is lower that the value found by Solanki et al.l ( l2003h . but the 
inversi ons favor the interpretation of the 3D loop reconstruc- 
tion of Solanki et al. (2003) over the interpretation proposed by 
iJudge (2009.) because the fitness of the results (a measure of the 
goodness of the fit) displays an asymmetric behavior, dropping 
rapidly towards lower heights, but comparatively gently towards 
greater heights. 

Further support for the above conclusion is provided by an 
inversion using HeLIx"*" of a profile at t he apex of eme rging 
loops in AR NOAA 10917 presented bv IXu et alj ( l2009l) . For 
that active region the inversion including atomic polarization 
returns a height of 7.1+2" (5.0 ± 1.4 Mm), which is n ot only 
higher than the formation height of 2.4 Mm favored bv iJudgd 
(l2009l) . but also higher than the reconstructed looptops at 4 Mm. 
Consequently, there is some discrepancy between the heights ob- 
tained from magnetic field reconstructions and from the Hanle 
effect, but in the two cases studied so far the difference has op- 
posite sign, which suggest that it is due to statistical errors. From 
the present analysis we cannot judge if the error is larger in 
the heights of the reconstructed loops or in the heights deduced 
from the Hanle effect. In both studied cases, however, the height 
deduced from the Hanle effect is higher than that proposed by 
I.Tudgel(l2009h . 

The interpretation put forward by I Judge! (l2009l) is probably 
valid in the quiet Sun, where under "normal conditions" the He 
1 10830 A multiplet tends to form in a thin layer at the top of 
the chromosphere. As Judge himself pointed out, the He triplet 
is formed over "a broad distribution of heights owing to the cor- 
rugated nature of the chromosphere". An example of a struc- 
ture in which He 1 10830 forms well above its quiet Sun forma- 
tion level are filaments or prominences. Thus, in prominences 
He 10830 A has been observed up to several tens of Mm above 
the limb. The reason is the much higher density of material at 
chromospheric temperatures inside prominences and filaments 
than in the surrounding corona. The density is sufficient to make 
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Fig. 4. The observed Stokes profile analyzed in this paper (solid 
black curves; left vertical scale; same reference direction as Fig. 
|2]i. The horizontal dashed line in the top panel corresponds to 
the adopted continuum level. Also plotted (red solid lines) are 
best fit profiles for a height of 6.3 Mm, obtained for a magnetic 
field of 320 G, with an inclination of 95° with respect to the local 
solar vertical, and an azimuth of -86° (see Section [33T l. 



the line optically thick. Our result suggests that the line forms 
at much greater heights than usual in an EFR, just as it does in 
filaments. Young, recently emerged loops in EFR, are often re- 
ferred to as arch filame nt systems ( AFS) when observed in Ha 
(see, e.g., the review bv lChoulll993 ) due to their arch-like shape 
and their high density of chromospheric material, which leads 
to enhanced absorption, giving them a filament-like appearance. 
The He triplet similarly displays enhanced absorption in both 
EFRs studied to date, displayin g dark filamen ts that look very 
similar to the ones seen in Ha IXu et alJl200"9l) . These dark fil- 
aments lie parallel to the magnetic field azimuth and are the lo- 
cations at which loops can be reconstructed. They do not show 
interruptions or inhomogeneities along their long axes as one 
would expect in the interpretation favoured by Judge, but rather a 
smooth change in brightness. The chromospheric magnetic field 
changes equally smoothly along these filaments, in stark con- 
trast to the photospheric field, which displays a very much higher 
level of inhomogeneity, including the pres ence of opposite po- 
larities b elow the l oop to ps (see Fig. 1 of ISolanki et all 120031 : 
Fig. 5 of lXuetal.Ll2009l) . Thus, the magnetic field structure in 
an EFR is rather different from the simple potential field model 
used by Judge (2009), in which the photospheric and chromo- 
spheric fields ar e basically identical. These points also favor the 
interpretation of lSolanki et alJ (l2003h and lend support to the re- 
constructed loops. 

Finally, in his paper I Judge! (12009!) did not put any empha- 
sis on the fact that the 3-D structure of the reconstructed loops 
agrees rather well with a non-linear force-free (NLFF) extrap- 
olation starting from a map of the photospheric vector mag- 
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netic field. This was derived from an independent inversion of Wiegelmann, T., Lagg, A., Solanki, S. K., Inhester. B., & Woch, J. 2005, A&A, 
simultaneously recorded spectropolarimetric measurements of 433, 701 

u . u • c- T A 1' c MI r-r- . 1 Xu, Z., Lagg, A,, & Solanki, S, K. 2009, A&A 

the photospheric Si I 10827 A line. Since NLFF extrapolations 
have no free parameter and take into account electric currents, 
this agreement is quite r emarkable. Interestingly the work of 
IWiegebnann et al.l (l2005h also showed that a potential field is 
a very poor representation of the chromospheric field, giving a 
completely different azimuth. Hence, even if ignoring the prob- 
lem that photospheric magnetic structure is far more complex 
than the chromospheric one, it is unclear if conclusions based on 
a (extremely simple) potential field computation (Judge 200^ 
have any validity for AR 945 1 . 

In summary, a state-of-art treatment of the Han le effect sug- 
gests that the interp retation of Solanki et al.l (l2003h is to be pre- 
ferred over that of I Judge! ( l2009h . More importantlv. this pa- 
per has demonstrated the value of treating the Hanle effect 
more comple t ely th an the simple approximation employed by 
ISolanki et alj (|2003|) . This more complete treatment leads not 
only to the possibility of estimating the height of formation of 
the He triplet, but also results in improved retrieved values of 
the magnetic vector. The effect of the height on the line profiles 
is rather small, however, and thus requires observations with ex- 
tremely low noise levels. 
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